ANL-HEP-CP-99-69 



X 



Open charm production in deep inelastic 
scattering at next-to-leading order at HERA * 

Brian W. Harris 

HEP Theory Group, Argonne National Laboratory, Argonne, IL 60439, USA 



Abstract. An introduction and overview of charm production in deep inelastic scat- 
tering at HERA is given. The existing next-to-leading order perturbative QCD calcu- 
lations are then reviewed, and key results are summarized. Finally, comparisons are 
Oh| made with the most recent HERA data, and unresolved issues are highlighted. 



1 Introduction 



^ . Electromagnetic interactions have long been used to study both hadronic struc- 

' ture and strong interaction dynamics. Examples include deep inelastic lepton- 

. nucleon scattering, hadroproduction of lepton pairs, the production of photons 

^\ ' with large transverse momenta, and various photoproduction processes involving 

. scattering of real or very low mass virtual photons from hadrons. In particular, 

heavy quark production in deep inelastic electron-proton scattering (DIS) is cal- 
culable in QCD and provides information on the gluonic content of the proton 
which is complementary to that obtained in direct photon production or struc- 
ture function scaling violation measurements. In addition, the scale of the hard 
^ i' scattering may be large relative to the mass of the charm quark, thus allowing 

O , one to study whether and when to treat the charm as a parton. 

Early measurements of open charm production in neutral current DIS, per- 
formed by the Berkeley-Fermilab-Princeton (BFP) and European Muon Col- 
laboration (EMC) experiments, touched upon the topics relevant for HERA 
today: production mechanism, charm fragmentation, gluon parton distribution 
' function extraction, and charm contribution to the proton structure function. 

See 11 for a review of these experiments. Activity in this area has increased 
recently with new data available from the HI Q and ZEUS j|] experiments at 
HERA. In particular, substantial samples of reconstructed D* hadrons have been 
obtained, and a semi- lepton decay mode analysis is underway]^. Considerably 
more data is anticipated in the next few years. 

Interest in the production mechanism is twofold. First, one is concerned with 
the leading twist-two term in the operator product expansion which incorporates 
the factorization theorem for hard scattering. Second, there is interest in studies 
of a higher twist charm component to the proton. 
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When considering the leading twist-two term the main issue is whether and 
when to treat the charm as a parton. Near threshold it is well established that 
charm is produced through photon-gluon fusion. On the other hand, very far 
above threshold the charm should be resummed into an effective parton dis- 
tribution. How one interpolates from one region to another is described by a 
variable flavor number scheme, several of which have been proposed recently 
must look at sufficiently inclusive observables in order to build up 
the logarithms that are to be summed, so predictive power for some observables 
is lost. In other words, for differential quantities, the fixed flavor number scheme 
of photon-gluon fusion provides the most appropriate formalism [jTo) . 

The idea of a higher twist charm component to the proton was introduced by 
Brodsky et al. . In this scenario the intrinsic charm quark Fock component in 
the proton wave function, \uudcc >, is generated by virtual interactions where 
gluons couple to two or more valence quarks. The probability for cc fluctuations 
to exist in a light hadron scales as 0(1/™^) relative to the twist-two component, 
and is therefore higher twist. The EMC data marginally support this idea. 
A full re-analysis of the EMC data has been carried out including both leading 
twist and intrinsic components at NLO including mass and threshold effects . 
The result is that one data point contributes to a 0.8 ± 0.6% normalization of 
the intrinsic component relative to leading twist. 

Fragmentation is the most contentious topic at present. The production and 
subsequent hadronization of charmed quarks in DIS is not as clean as the much 
studied case of production in e^e~ annihilation. In particular, in DIS there are 
proton beam remnants around which necessarily talk to the charm quark as it 
hadronizes. If events too close to the beam direction are selected, one can ex- 
pect deviations from models which do not account for this effect. Presently it is 
hoped that detailed experimental studies in the Breit frame, wherein one hemi- 
sphere resembles e+e^, will provide additional information on the fragmentation 
process. 

Because the production method is dominated by the leading twist photon- 
gluon fusion near threshold, it is possible to extract a gluon parton distribution 
function (PDF) and compare with existing PDFs, which derive their gluon in- 
formation from other sources, such as direct photon production or structure 
function scaling violation measurements. This has been done by both HI and 
ZEUS collaborations. The results are completely consistent with comparable 
PDF sets in the same x and range. While unlikely to ever replace structure 
function scaling violation measurements, charm production does serve as a nice 
consistency check on the gluon PDF. 

One may also wish to take the open charm measurement and extrapolate the 
cross section over the full phase space and then extract the charm's contribution 
to the proton structure function. Historically the structure functions have been 
useful as input into global analyses and a testing ground for the variable flavor 
number schemes mentioned above. Care must be taken so that conclusions drawn 
are not artifacts of the extrapolation procedure (model). 
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Fig. 1. Typical Feynman diagrams contributing to the amplitude for neutral current 
charm production, (a) An order egs Born diagram, (b) An order egl virtual diagram, 
(c) An order egl gluon-bremsstrahlung diagram, (d) An order eg^ light quark initiated 
diagram 



2 Next-to-leading Order Calculations 



The reaction under consideration is charm quark production via neutral-current 
electron-proton scattering, {l)+P{p) {l') + c{pi)+X . When the momen- 

tum transfer squared = —q^ > (g = / — /') is not too large <^ the 
contribution from Z boson exchange is suppressed and the process is dominated 
by virtual-photon exchange. After an azimuthal integration, the cross section 



may be written in terms of structure functions i^2^(a;, Q 
as follows: 



and Fl{x, Q , rric 



{ [1 + (1 - yf] F^{x, Q\ m,) - y^Fl{x, Q^ m,)} (1) 



where x = /2p-q and y — p-q/p-l are the usual Bjorken scaling variables, a is 
the electromagnetic coupling, and rric is the charm mass. The scaling variables 
are related to the square of the center of momentum energy of the electron- 
proton system S — {I + p)'^ via xyS — Q^. The total cross section is given by 

m 



dy 



dQ^ 



\dydQ^ 



(2) 



where me is the electron mass. 

Typical Feynman diagrams for this process are shown in Fig. ^ The inter- 
ference of the 0{egs) Born diagrams with the ©(egf ) one-loop virtual diagrams 
produces a result that is 0{aa^). The ultraviolet divergences are removed by 
renormalization in the Collins- Wilczek-Zee scheme The result is added to 
the square of the gluon-bremsstrahlung diagrams which is also 0{aa^). Initial 
state coUinear singularities are mass factorized to obtain a finite NLO result. 
For the light quark initiated subprocess there are no virtual contributions at 
0{aaj.). One only encounters initial state coUinear singularities which are again 
removed using factorization. The full NLO corrections were first calculated in 
[fsl and may be written in the form 
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l=q.,q 

,2 



(3) 



with k — 2, L. The lower limit on the integration over the parton momentum 
fraction ^ is ^min — x(4m^ + Q^)/Q^- The parton momentum distributions in 
the proton are denoted by fi/p{^, ^J,'^)- The sum is taken over the light quarks, 
q = u,d, s. The mass factorization scale /i/ has been set equal to the renormal- 
ization scale /i^ and is denoted by fj,. All charges are in units of e. cfl, c^'^], 

cl^l , (i — g,q,q), and Sf^l, o\^l , (i = q,q) are scale independent parton coeffi- 
cient functions, and are distinguished by their origin. The c-coefficient functions 
originate from processes involving the virtual photon-heavy quark coupling, the 
d-coefficient functions arise from processes involving the virtual photon-light 
quark coupling, and the o-coefficient functions are from the interference between 
these processes. 

In addition to describing the structure functions, the coefficient functions 
of (l^ may also be used to compute the single-inclusive distributions dFk/dpt 
and dFk/dy [ p^ , where pt and y are the transverse momentum and rapidity, re- 
spectively, of the heavy quark in the virtual photon-proton center of momentum 
system. 

Unfortunately, the analytic expressions of the coefficient functions are too 
long to be published in journal form. They are however available as computer 
code, but, for the same reason, the code tends to be slow. Initially this was seen 
as an impediment to including them in a global fitting program. Then in ||l7|| , 
with the threshold and asymptotic behavior of the coefficient functions removed, 
it was possible to numerically tabulate grids, with a fast interpolation routine, 
so that speedy computation of Eq. - (||) became possible. This process was 
further refined in [T^ . 

The coefficient functions in the large momentum transfer limit are a necessary 
ingredient for constructing the variable flavor number schemes mentioned in the 
introduction. Exact analytic formulae for the d^,^] together with analytic formulae 
for the coeflacient functions cf], c\^l, c\^l ,(i = g, q, q), and d'j^] ,{i — q, q) in the 
limit ^ can be found in Q. 

The coefficient functions have also been calculated in a fully differential form 
[ pO| . These in turn can be used with Eq. (^) to construct pair- inclusive distri- 
butions such as dFk / dMcs (2|] where M^g is the invariant mass of the produced 
charm-anticharm system. 
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The resulting differential structure functions and Eq. have further been 
used to construct a NLO monte carlo style program, HVQDIS |22,2^. The 
basic components (in terms of virtual-photon-proton scattering) are the 2 to 2 
body squared matrix elements through one-loop order and tree level 2 to 3 body 
squared matrix elements, for both photon-gluon and photon-light-quark initiated 
subprocesses, as shown in Fig. |l|. It is therefore possible to study fully-, single-, 
and semi-inclusive production at NLO, and three body final states at leading 
order. The goal of this NLO calculation is to organize the soft and coUinear 
singularity cancellations without loss of information in terms of observables that 
can be predicted. 

The subtraction method provides a mechanism for this cancellation. It allows 
one to isolate the soft and coUinear singularities of the 2 to 3 body processes 
within the framework of dimensional regularization without calculating all the 
phase space integrals in a space-time dimension n ^ A. Expressions for the three- 
body squared matrix elements in the limit where an emitted gluon is soft appear 
in a factorized form where poles e = 2 — n/2 multiply leading order squared 
matrix elements. These soft singularities cancel upon addition of the interference 
of the leading order diagrams with the renormalized one-loop virtual diagrams. 
The remaining singularities are initial state coUinear in origin. The three-body 
squared matrix elements appear in a factorized form, where poles in e multiply 
splitting functions convolved with leading order squared matrix elements. These 
coUinear singularities are removed through mass factorization. 

The result of this calculation is an expression that is finite in four-dimensional 
space time. One can compute all phase space integrations using standard monte 
carlo integration techniques. The final result is a program which returns parton 
kinematic configurations and their corresponding weights, accurate to 0{aa1). 
The user is free to histogram any set of infrared-safe observables and apply cuts, 
all in a single histogramming subroutine. Additionally, one may study heavy 
hadrons using the Peterson et al. model |2j]. Detailed physics results from this 
program and a description of the necessary cross checks the program satisfies 
are given in [E3|. See also for recent improvements. 



3 Results and Current Issues 

Charmed meson differential cross sections are measured experimentally 
within some detector acceptance region, and the corresponding theory predic- 
tions can be made using HVQDIS. As we saw above, the cross section is an 
integral over the structure functions. Therefore, the two share many of the same 
properties. Before discussing the cross sections, the more salient features of the 
NLO structure functions will be reviewed. The interested reader can find addi- 
tional details in the original paper [ psf and, more so, in the recent phenomeno- 
logical analyses [|8U26H^,|||. 



For moderate ~ 10 GeV^ one finds that the charm quark contribution at 
small X ^ 10~^ is approximately 25% of the total structure function (defined as 
light parton plus heavy quark contributions). In contrast, the contribution from 
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Fig. 2. The renormalization/factorization scale dependence of the structure function 
FS{x,Q^,mc) for = 3 GeV^ (left) and = 25 GeV^ (right) for various x values. 
The results for GRV94 {solid lines) and GTEQ4F3 {dashed lines) parton distribution 
sets are shown 



bottom quarks is only a few percent due to charge and phase space suppression. 
The structure functions show a marked rise at small x due primarily to the 
rapidly rising gluon distribution: the gluon initiated contributions comprise most 
of the structure function. The light quark initiated processes give only a few 
percent contribution at small x. The scale dependence of the structure functions 
is very small in the HERA x and regions. This is demonstrated in Fig. || 
for various x and values. The largest variation comes from our imprecise 
knowledge of the charm quark mass. For example, a ±10% variation of the 
charm mass about the central value of 1.5 GeV gives a ±20% variation in the 
structure function for small x and moderate Q^. 

At moderate Q^, and x values larger than 0.01, the charm structure func- 
tion is increasingly dominated by partonic processes near the charm quark pair- 
production threshold. The large size of the gluon density for small momentum 



fractions gives relatively large weight to such processes |27|. Although the QCD 
corrections at presently accessible x values are moderate (about 30 — 40%), with 
an increasing amount of data to be gathered at higher x, it is worthwhile to have 
a closer look at such near-threshold subprocesses. In this kinematic region, the 
QCD corrections are dominated by large Sudakov double logarithms. Recently 
[ p9| , these Sudakov logarithms have been resummed to all orders of perturbation 
theory, to next-to-leading logarithmic accuracy, and, moreover, in single-particle 
inclusive kinematics ||30| . Let us recall the main results. First, the quality of the 
approximation for the neii-to-leading logarithmic threshold resummation was 
found to be clearly superior to leading logarithmic one. Furthermore, the resum- 
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Fig. 3. Next-to-leading order diflterential cross sections for charm production at ^/S = 
301 GeV at HERA using the GRV94 {dashed lines) and CTEQ4F3 [solid lines) parton 
distribution sets, (a) Transverse momentum, (b) Pseudo-rapidity, (c) Bjorken x {left 
set) and parton momentum fraction ^ [right set), (d) Photon virtuality 



mation provided next-to- next-to-leading order estimates [g9| , which were found 
to be sizable for x > 0.05. 

Let us now examine some of the properties of the charm quark cross section. 
Results are presented in the HERA lab frame with positive rapidity in the proton 
direction. The proton and electron beam energies are taken to be 820 GeV and 
27.6 GeV, respectively. Results are for the kinematic range 3 < < 50 GeV^ 
and 0.1 < y < 0.7. The CTEQ4F3 [|l) and GRV94 HO || proton-parton 
distribution sets are used. The renormalization and factorization scales have 
been set equal to ^. 

Fig. H shows the NLO cross sections differential in transverse momentum pt, 
pseudo-rapidity 77, Bjorken x, and momentum transfer for charm quark pro- 
duction using the GRV94 (dashed) and CTEQ4F3 (sohd) parton distribution 
sets at fj, = a/Q2 _|_ 4jy^2 j^i^ — 15 GeV. From Eq. ^ the parton distri- 

butions are probed at a momentum fraction ^ which is typically one order of 
magnitude larger the x. This is illustrated in Fig. ^ where a plot of dcT/dlog(^) 
vs. log(^) (right set of curves) is superimposed on the plot of da/d\og{x) vs. 
log(a;) (left set of curves). The difference between the curves is approximately 
10% at ^ = 10-2-^. 

The scale dependence of the NLO differential cross sections is shown in Fig. ^ 
The curves were made using the CTEQ4F3 parton distribution set at /i = 2mc 
(solid) and /i = 2^/Q^ + 4m^ (dashed) with nic = 1.5 GeV. The curves show 
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Fig. 4. Same set of distributions as Fig. ^ but this time showing the variation 
with respect to renormalization/factorization scale, ^ = 2mc {solid lines) and jj, = 
2yQ^+4mf {dashed lines) 



very little scale dependence. This can be anticipated from the results shown in 
Fig. 1^ and the distribution in Bjorken x shown in Fig. The latter shows the 
cross section is dominated by x ^ lO^'^'^ — 3.2 x 10~^ while the former shows 
that, independent of Q^, the structure function is very flat in this particular x 
region. Therefore, the cross section tends to be fairly insensitive to the choice of 
scale. Other kinematic regions show increased scale dependence. 

The largest uncertainty in the structure function calculation is due to the 
charm quark mass. The same is true for the cross section as shown in Fig. ^. The 
NLO differential cross sections for charm quark production using the CTEQ4F3 
parton distribution set at fi — \/Q'^ + Arn^ with nic ~ 1.35 GeV (solid) and 
rric = 1.65 GeV (dash) are shown. Mass effects are smaller at the larger transverse 
mass because they are suppressed by powers of rric/pt in the matrix elements. 
However, if the range is extended much further, large logarithms of the form 
ln(p( /rn^) appear in the cross section and should be resummed. 

Before closing, we compare the NLO calculations described above with the 
most recent HERA data from HI Q] and ZEUS |^ collaborations. The mea- 
surements make use of a tagging technique wherein the D* meson kinematics 
are reconstructed using the tight constrains of the decay D*^ -^^''^tiow ~^ 
{K^'K^)n'^i^^. In order to make the comparison, the theory prediction for the 
charm production cross sections must be converted to those of charmed meson 
production. This is done using a simple non-evolving Peterson et al. model 
which depends on one parameter e which is taken from e+e^ data. The overall 
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Fig. 5. Same set of distributions as Fig. |3|, but this time showing the variation with 
respect to charm quark mass, rric = 1.35 GeV {solid lines) and rric = 1.65 GeV {dashed 
lines ) 



cross section normalization is set by the hadronization fraction /(c —>■ D*) again 
taken from e^e" data. The four-vector for the D* is constructed from that of 
the charmed quark by smearing the the charm three-vector in the lab frame. 
The energy component is then fixed such that the four-vector has the physical 
D* mass, 2.01 GeV. 

Shown in Fig. || are D* meson cross sections measured by the HI collabora- 
tion Q differential in transverse momentum pseudo-rapidity ?], reconstructed 
parton momentum fraction x°^^ ^ and momentum transfer compared to the 
NLO calculation described in the Sec. ^ plus the Peterson et al. fragmentation 
model. The shaded band corresponds to varying the charm quark mass of 1.5 
GeV by ±0.2 GeV. Overall the agreement is good, except for the pseudo-rapidity 
plot in which the theory under (over) estimates the data in the forward (back- 
ward) region. 

The ZEUS collaboration ||^ has also measured D* meson cross sections 
differential in momentum transfer Q^, Bjorken x, hadronic energy W, trans- 
verse momentum pt, pseudo-rapidity 77, and D* momentum fraction x{D*) = 
2|P7Pcins|/W^ which are compared with theory in Fig. ^ The boundaries of the 
bands correspond to varying the charm quark mass of 1.4 GeV by ±0.1 GeV. 
Again, the overall agreement is good, but the theory underestimates the data 
in the forward region and overestimates it in the backward region. Additionally, 
the D* momentum fraction data, which is particularly sensitive to the charm 
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Fig. 6. Various differential cross sections for D* meson production at HERA as mea- 
sured by the HI collaboration |^] compared to the next-to- leading order calculation 
described in the Sec. |^ plus a Peterson et al. fragmentation model. The shaded band 
corresponds to varying the charm quark mass from 1.3 to 1.7 GeV 



hadronization process, is poorly described. Similar effects are seen in the D* 
photoproduction data at HERA. 

Variations of the parton distribution set, renormalization/factorization scale, 
charm mass, or fragmentation parameter e are unable to account for the differ- 
ences between data and theory. It also appears unlikely that an evolving frag- 
mentation function would help; and the momentum transfers are large enough 
that any photon structure is surely negligible. 

One explanation |]33| proposed for the photoproduction data appears to work 
for the DIS data as well. Qualitatively, one is invited to think of a color string 
connecting the hadronizing charm quark and the proton remnant which pulls 
(drags) the charmed meson to the forward region. This is made quantitative in 
the Lund String model modified for heavy flavor production , as implemented 
in Pythia |^ . The shaded band in Fig. |^ shows what happens when the Peterson 
et al. model is replaced by an effective fragmentation model extracted from the 
Pythia based monte carlo RAPGAP|^. The agreement is much improved. 

Another way to improve the agreement between data and theory is to simply 
raise the minimum of the events that are selected. Data from a slightly different 
decay chain, but higher minimum pt cut are shown as open triangles. Here the 
Peterson and RAPGAP improved NLO predictions give essentially the same 
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Fig. 7. Various differential cross sections for D* meson production at HERA as mea- 
sured by the ZEUS collaboration compared to the next-to-leading order calculation 
described in the Sec. ^ plus a Peterson et al. fragmentation model. For the shaded 
band the Peterson et al. model was replaced by an effective fragmentation model ex- 
tracted from RAPGAP which includes a drag effect between the proton remnant and 
the hadronizing charm. The bands result from varying the charm quark mass from 1.3 
to 1.5 GeV 



results, as expected. However, a seemingly large fluctuation in the forward most 
data bin somewhat clouds this expectation. 

A number of additional studies have been done. For example, HI us- 
ing the above cross sections, has extracted a NLO gluon PDF which agrees 
well with their own gluon PDF, obtained indirectly through structure function 
scaling violations, and that of CTEQ4F3. In a different approach, ZEUS 
has extrapolated over the full phase space and extracted the structure function 

^charm 

In closing, the next-to-leading order calculations described herein have been 
very successful in describing charm production at HERA. A variety of different 
observables have been studied, and the gluon parton distribution function and 
the charm contribution to proton structure function have been extracted. The 
weakest stage of the calculation is, not surprisingly, modeling the hadronization 
of the produced charm to the observed charmed meson, especially at low trans- 
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verse momentum. One could take this as an opportunity to study hadronization 
in the presence of beam remnants. 
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